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Abstract 

^' 
Ch ' The ground-state phase diagram of an Ising spin-glass model on a random graph 

^ , with an arbitrary fraction (w) of ferromagnetic interactions is analysed in the pres- 

ence of an external field. Using the replica method, and performing an analysis of 
stability of the replica-symmetric solution, it was shown that w = 1/2, correponding 

^ . to an unbiased spin glass, is a singular point in the phase diagram, separating a re- 

'nT I gion with a spin glass phase {w < 1/2) from a region with spin glass, ferromagnetic, 

mixed, and paramagnetic phases {w > 1/2). 

O ■ PACS numbers: 75.10.Nr, 75.50.Lk, 89.75.-k 

O 

The mean-field formulations of the Ising spin glass and related disordered models on 

a lattice are known to lead to replica - symmetry breaking, which is in turn associated 

with the existence of a large number of metastable states. The number of local minima 

'^ I in the free energy landscape of these solutions grows exponentially with the size of the 

Q I lattice^] . The impressive richness of these mean-field results has been a motivation to look 

O ' at more realistic spin-glass models, with the inclusion of the effects of either short-range 

interactions or the finite coordination of the crystal lattices. Besides the implementation of 

K> ! the Ising spin glass (SG) on a Bethe lattice j2], which is one of the natural ways to deal with 

H I a finite connectivity, there is a work by Viana and BrayjH] for the Ising SG on a random 

graph with a mean connectivity. According to a number of investigations|ll 00111, this 
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model of Viana and Bray (VB), which is tractable by the replica method, leads to a glassy 
phase in the ground state. However, due to the subtleties of the problem, there are just a 
few available analytical results, and a complete solution remains an open question. Also, 
in spite of some efforts, there are no satisfactory solutions for the field behavior of the 
SG problem on random graphsjS]. In contrast to earlier results, more recent analytical 
calculations for the ground state of a dilute ferromagnet, with a specific choice of random 
fields, do not indicate the existence of replica breaking P, TU]. 

In this article, some results for the ground state of a frustrated ± J Ising spin-glass on a 
random graph, in the presence of an external magnetic field, and with an arbitrary fraction 
of ferromagnetic interactions are reported. A parameter w that gauges the concentration 
of these ferromagnetic (Jq > 0) interactions is introduced. For w = 1/2, this model is 
a pure, unbiased, spin glass; ferromagnetic (antiferromagnetic) interactions are favoured 
for w > 1/2 {w < 1/2). It is known that the presence of an external field always increases 
the technical complexity of the solutions of these problems. In this analysis, the phase 
diagram in the ground state for external fields that are multiples of Jo is obtained, which 
resembles the choice of random fields that was already adopted in a previous work jH] . The 
topology of the random network favours the appearance of a spin glass phase for w < 1/2. 
For w > 1/2, the usual spin glass, ferromagnetic, mixed, and paramagnetic phases are 
obtained. It should be emphasized that even in the presence of overall antiferromagnetic 
interactions {w < 1/2) the inherent disorder of the lattice still yields a spin-glass phase. 

Consider the Hamiltonian 

T^iWx}) = - 2_^ Jx,x'CrxO'x' - H 2_^ CTa;, (1) 

(a;,a;')gAjv xSAjv 

where a^ G { — 1,1} is an Ising spin on the site x, belonging to a finite lattice A^r of A^ 
sites and H is the external magnetic field. The first sum is over all distinct pairs of lattice 
sites. The exchange integrals {Jx,x'} are independent and identically distributed random 
variables, associated with the distribution 

Pj{Jx,x') = y-~ J^J ^iJx,x') + J^p{Jx,x') , (2) 

where 

p{Jx,x') = wS{Jx,x' - Jo) + (1 - w)5{Jx,x' + Jo) , (3) 

and Jo is a positive parameter. According to the distribution Q, two sites are connected 
(disconnected) with a probability c/N (1 — c/N), where c is the mean connectivity of 
the lattice. In this diluted ± J model, parameter w weighs the fraction of ferromagnetic 
bonds. 



Using the replica method, the variational free energy is written as 

if ^ 

^ (^ r=0 (ai,...,a,) 

n 

G'(K})+/3/7 5^a, 



In Trjo-} exp 



a=l 



(4) 



where 



h, = / dJp{J) cosh"(/3J) tanh''(/3J) , (5) 



and 

n 

GiWa}) = C^ J^ &rgai,-,a,Cr«^ ■ ■ ■ f^a, (6) 

r=0 (ai,--- ,ar) 

is the global order parameter [^ ITT] . which takes into account the 2" order parameters of 
this problem (in contrast to the Sherrington - Kirkpatrick (SK) modeipp, which has just 
two order parameters). 

The stationary free energy comes from the minimization of the variational free energy 
(jlj) with respect to the set {qai,-,ar}y which leads to the stationary conditions 

^ Tri^l^T^^ ■ ■ ■ g^,. exp [^({0-^}) + f3H J2l^-^ a^] 
^""■■■'"^ Tr.exp[G(K})+/3/7ELi^J ' 

The replica-symmetric (RS) solution of this problem may be obtained from the intro- 
duction of a local field h, with an effective probability distribution P{h), which is equally 
applied to all of the replica spin variables. Thus, one has 



Qai,---,ar= dhP{h) tanlf (Ph) . 
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Moreover, in the replica - symmetric context, one can see that G{{cx}) = G{a), where 

Equation (jH} leads to a relation between the distribution of local fields, P, and the 
global order parameter, G. The analysis will be carried out at the ground state. Also, 



in analogy to a previous workjH], the RS solutions of this problem are restricted to the 
condition r = H/ Jq, where r is an integer. In other words, the solutions in this work are 
written in the form 



P{h) = Y,cik6{h-kJ^), 



(9) 



kei 



which takes into account the discrete part of the distribution function only; the coefficients 
{a/c} are given below. Under these conditions, one can show that 



where 



G{iy) = lim G{iy/p) =A + B'e'y^° + C'e"'^^" 



A = cao , B' = wB+{l-w)C , C = {I - w)B + wC , 



(10) 
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with 



B = c y^ ttfc , C = c y^ ttk and ak = e 



A-c 



k=l 



k=—oc 



k — r 

B'\~ 

'C' 



Ik-r{2VB'C' 



'121 



where Iu{x) is the modified Bessel function of order z/. 

The solution of this problem depends on three equations, from which one can obtain 
A, B and C. Using the previous equations, one finds 



C 



B' 



A = c[ — ) Iri2^/WC') , c = A + B' + C'{ = A + B + C) 



(13) 



and 



B' 



ce 



A-c 



, — 4 r" r- 

we - we Oro- 



r-l 1 

-wC"(^J I due^'"" (1 - u)"^ /,_! {2^B'C'{1 - u) 

Ir (2^/B'C'{l-u 



+ {l-w)C'(^X I due''''' a- u)-^ 



(14) 



The stability analysis can be performed along the same steps of previous works jHl El- 
It is based on the search of the eigenvalues of the Hessian matrix associated with the 
variational free energy. By some technical manipulations of the eigenvalue equations, 
which are similar to previous calculations |1H IH]. a separate analysis for the cases w ^ 1/2 
and If = 1/2 is carried out. 

For w 7^ 1/2, one finds the following transversal eigenvalues: 

I A I A 1 , A \ (w , 

and — + A 



c ' c c ' c c(l — 2w) c I — 2w \ c 

(15) 



where A± = ±J [^ ± ai (1 - 2w)] [^ ± a_i (1 - 2w)] . The stability of the RS phase 
depends on the smallest eigenvalue, 

. _ A w 1 ( A- , w>l ,^„. 

^'-~^~ c{l-2w)'l^2^''\ A+ , w<\ ■ ^^^^ 

Some limiting cases are analytically accessible. In the pure ferromagnetic case {w = 1), 
one has \q{w = 1) = — ^ + c ~ v^QiO-i, which confirms the stability of the RS solution in 
large external fields, in which case Ao(w = 1, if ~ oo) ~ 1/c (since, for a fixed value of fc, 
Qk goes to zero for sufficiently large H). Moreover, by the same argument, one can make 
Ao positive for any w > 1/2 if the field is sufficiently large. On the other hand, the w = 



(pure antiferromagnetic interactions) case implies Xo{w = 0) = — 7 — y [^ + o-i] [^ + ct-i] , 
which is clearly negative. This means that the RS solution for this finitely connected 
spin-glass model is always unstable in this regime, even in the limit H — > 00, which yields 
Xo{w = 0,H ^ 00) = —1/c. Note that in the limit of infinite connectivity, c -^ 00, this 
eigenvalue goes to zero, becoming marginally stable as in the SK model. 

The phase diagrams (Figure 1), which are qualitatively similar for other choices of the 
mean connectivity c, show that the RS solution is always unstable for w < 1/2, which 
can also be seen from equation (fTBj) . For w > 1/2, the RS solution is stable above a 
critical value w* of the concentration, which depends on the applied field. This behaviour 
can be illustrated in the r x w phase diagram, where the SG phase is separated from the 
paramagnetic phase by a de Almeida-Thouless line|12j. 

In zero field one can show the presence of ferromagnetic and mixed phases. This mixed 
phase is characterized by the asymmetric distribution of the local field, which leads to 
ttk 7^ a_fc for some values of A; G Z. 
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Figure 1: Left: Graphic Xq x w, obtained by numerical calculation of the eigenvalue. 
Right: Spin glass (SG), paramagnetic (P), ferromagnetic (F) and mixed (M) phases in 
r X w phase diagram. 



For w = 1/2, the eigenvalues of the stability matrix are 
1 



A 



A 1 , All, 

h - and h - ± - (cti + ct-ij 

c c c c 2 



(17) 



The smallest of these eigenvalues is given by 



An 



A 1 

'- + - 

c c 



(cti + Ct- 



This eigenvalue does not have a definite sign, but it is positive (tending to 1/c) for 
sufficiently large field, and agrees with equation (fTHj) when w I 1/2 (Ao(w) is continuous 
from right). On the other hand, it is easy to see that lim^|i/2 Xo{w) = — oo. 

In summary, this work reported the ground-state phase diagram of a diluted ± J Ising 
SG model, in an external field, with a fraction w of ferromagnetic bonds. For w < 1/2, 
which corresponds to a larger concentration of antiferromagnetic bonds, the topology of 
the lattice leads to the existence of a spin-glass phase only. This behaviour has also been 
found numerically for free scale networks with antiferromagnetic interactions only [THj . 
For w > 1/2, with a larger concentration of ferromagnetic bonds, the usual phases are 
found. 
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